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Abstract 

We calculate rapidity and transverse momentum distributions of charmonium formed in high 
energy heavy ion collisions from incoherent recombination of charm quarks. The results are very 
sensitive to the corresponding distributions of the charm quarks, and thus can serve as a probe of 
the state of matter produced in the heavy ion collision. At one extreme we generate a set of charm 
pair momenta directly from pQCD amplitudes, which are appropriate if one can neglect interaction 
of the quarks with the medium. At the other extreme we generate momenta of charm quarks in 
thermal equilibrium with the expanding medium, appropriate for an extremely strong interaction. 
Explicit predictions are made for J/i/j formation in Au-Au interactions at RHIC. We find that for 
the case in which charm quark momenta are unchanged from the pQCD production calculation, 
both the rapidity and transverse momentum spectra of the formed J/ip are substantially narrower 
than would be anticipated in scenarios which do not include the in-medium formation. In partic- 
ular, the average transverse momentum of the J ftp will exhibit a non-monotonic behavior in the 
progression from p-p to p-A to A- A interactions. 
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I. INTRODUCTION 



The behavior of heavy quarkonium states as an experimental signature of medium modi- 
: ication in high energy heavy ion collisions originated with the prediction of J / ip suppression 
lj], which follows from screening of the color-confinement potential above the phase transi- 
tion temperature. The suppression argument requires that the probability for recombination 
of the c and c quark is negligible. Recently it has been realized that for RHIC and LHC 
collider energies, there can be a modification of the suppression argument DD 

Consider a central heavy ion collision in which multiple cc pairs are produced in binary 
nucleon-nucleon interactions. If these quarks are then immersed in a medium (such as a 
quark-gluon plasma) in which a charm quark from one initial cc pair can readily interact 
with an anticharm quark from a different initial cc pair, one expects that the number of J / ip 
formed from such interactions will be proportional to the total number of cc combinations. 
The relative efficiency of this quadratic process must be normalized by the number of inter- 
actions involving charm plus light-flavor quarks in which open-charm hadrons are formed. 
The normalization measure is generally taken [4] to be proportional to the total number of 
charged hadrons produced in the heavy ion collision, N ch , which leads to Nj/^ oc N c5 2 /N ch . 
(Of course this behavior must saturate before Nj/^ becomes comparable in magnitude to 
N c5 .) Then at sufficiently large N c5 this quadratic behavior must dominate over the color- 
screening suppression. 

Initial estimates pj from extrapolations of the low energy charm production measure- 
ments and predictions of pQCD indicated of N c5 « 10 for central RHIC collisions and 
several hundred at LHC. Measurements at RHIC of high transverse momentum leptons in 

mesons in d-Au collisions by STAR [7J require N cS ~ 40. Statistical and systematic uncer- 
tainties are large enough such that these two measurements are consistent, but the general 
trend indicates a more rapid growth with energy than initially estimated [q]. 

Two distinct realizations of this mechanism for heavy quarkonium formation have been 
developed. In the statistical hadronization model, quarks are distributed into hadrons dur- 
ing the hadronization transition according to chemical equilibrium ratios. The total number 
of heavy quarks is assumed to be determined by initial production, and an extra fugacity 
factor 7/ is determined by overall flavor conservation. In the case of hidden flavor such as 
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J/if>, two powers of 7 C modify the thermal equilibrium density, which leads to the expected 



quadratic dependence on N c5 |9 



101 ] . There is no explicit prediction for the J/ijj momentum 



spectrum, but it is natural to assume a thermal distribution in this model at the decon- 
fmement temperature T c « 170 MeV. It is noteworthy that ratios involving radially-excited 
states such as -rrr retain their unmodified statistical values. 

n n 



The kinetic formation model 3j, |4 , 
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12, 



13J assumes that heavy quarkonium formation 



takes place during the entire lifetime of a color-deconfmed phase. Predictions for the result- 
ing population require the specification of cross sections for formation and dissociation in 
the medium. The motivation for such a scenario in the case of J /if) formation has received 
support from recent lattice calculations of spectral functions. These indicate that the J /if) 
and r] c will exist in a thermal environment at temperatures well above the deconfinement 



transition 



14. 



Note that this property is likely valid only for the most deeply bound 
states. Thus, unlike the statistical model, one does not expect a simple relation to exist 
between the J /if) and if)' population. Also in contrast with statistical hadronization, the 
momentum spectrum of formed J /if} will reflect the initially-produced charm quark spectra 
plus any modification due to interaction with medium. 

In principle, both of these formation processes could occur sequentially. Considerable 
progress has been made H Q 

in such a scenario when both bound and unbound 
charm populations evolve in a thermal fireball. 

In the next section we review the basic properties of kinetic in-medium formation of heavy 
quarkonium, using a specific model for J/ if) which exhibits the expected general properties. 
We have previously found that the initial PHENIX data from RHIC is able to constrain 
the model parameters within a fairly broad range, but is unable to confirm or rule out 
this type of formation. A more detailed test must await comparison with the momentum 
spectra of the formed J/if). Section ITTTl contains details of the two specific charm quark 
distributions which we consider. One is calculated directly from the initial production 
process using NLO pQCD amplitudes to generate a set of cc pair events appropriate for 
RHIC energy. The other generates a corresponding set of events which would be expected 
if the charm quarks subsequently come into local thermal equilibrium with the expanding 
medium before forming the J /if). In Section HVl we present a detailed model calculation of the 
normalized rapidity and transverse momentum spectra for in-medium J/ if) formation, using 
charm quark spectra calculated from the pQCD amplitudes. It is important to keep in mind 
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that we are searching for a signature in the J/ ip rapidity and transverse momentum spectra 
which would indicate the presence of in-medium formation, and which is independent of the 
absolute magnitude of the formation process. Our results are shown to be quite robust with 
respect to substantial variations in the model parameters. Section IV Al contains a specific 
for RHIC results, where we use data for pp and dAu interactions to constrain some of the 
model parameters. Section IV Bl presents the corresponding calculations using charm quarks 
in local thermal equilibrium. Comparisons with predictions of other models and a summary 
of results complete this presentation. 



II. FORMATION OF J/ijj IN A COLOR-DECONFINED MEDIUM 

We first examine the net number of J/ip produced in a color-deconfined medium due 
to the competing reactions of (a) formation involving recombination of c and c and (b) 
dissociation of J/ip induced by interactions with constituents of the medium. The simplest 
dissociation reaction utilizes absorption of single deconfined gluons in the medium to ionize 
the color singlet J/ift, g + J/ift—tc+c, resulting in a cc pair in a color octet state. This 
process was originally proposed 20| as a dynamical counterpart of the static color screening 
effect. The inverse of this process then serves as the corresponding formation reaction, in 
which a cc pair in a color octet state emits a color octet gluon and falls into the color singlet 
J/ijj bound state. 

One can then follow the time evolution of charm quark and charmonium numbers accord- 
ing to a Boltzmann equation in which the formation and dissociation reactions compete. 

= A F iV c N.iVit)]- 1 - X D N J/lP Pg , (1) 

with p g the number density of gluons in the medium. The reactivity A is the product of 
the reaction cross section and initial relative velocity (av re \) averaged over the momentum 
distribution of the initial participants, i.e. c and c for Xp and J/ip and g for Xp. The gluon 
density is determined by the equilibrium value in the medium at each temperature, and the 
volume must be modeled according to the expansion and cooling profiles of the heavy ion 
interaction region. 

This equation has an analytic solution in the case where the total number of formed J/ ip 
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is much smaller than the initial number of N r p- 



f 



Nj/^tf) = e{t f )[N m {t ) + N% / A F [V(t) e(t)]- 1 dt], (2) 

J t 

where to an d tf define the lifetime of the deconfined region. Note that the function e(t/) = 
e ~Jt r>p g dt wou \d De the suppression factor in this scenario if the formation mechanism 
were neglected. 

The initial calculation, Q used the r atio of nucleon participants to participant density 
calculated in a Glauber model to define a transverse area enclosed by the boundary of the 
region of color deconfinement. This is supplemented by longitudinal expansion starting at 
an initial time t Q = 0.5 fm. Transverse expansion was initially neglected, but has been 
included in subsequent calculations The expansion was taken to be isentropic, which 
determines the time evolution behavior of the temperature T(t). The initial value T is 
taken as a parameter, and the final Tf is fixed at the hadronization point. 

The reactivities A^ and A^ require specification of cross sections. For o D we use the 
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22, 



OPE-based model of gluon dissociation of deeply-bound heavy quarkonium 
which is related via detailed balance to the corresponding a p. Written in terms of the heavy 
quarkonium reduced mass //(= mg/2), gluon momentum in the cc rest frame k and binding 
energy e , the dissociation cross section is 

_ 2tt /32\ 2 / 2/ A 1/2 1 (fc/ep - l) 3 / 2 

This expression assumes that the heavy quarkonium has a nonrelativistic Coulomb bound 
state spectrum with eo >> Aqcd, and utilizes an operator product expansion in the large 
A^ c limit. These cross sections for J/i(j kinematics are shown in Fig. ^ One sees that they 
are peaked at low energy, and that op > <Jp> due to the large binding energy (we used the 
vacuum value e ~ 600 MeV). 

We should note that the approximations concerning a nonrelativistic Coulomb bound 
state spectrum are somewhat marginal for the charmonium case, although they are better- 
justified for the upsilon states. Therefore for applications to J/ip we will investigate the 
effects of a considerable range of alternative cross section behavior, as has been done in a 



number of treatments of J /if} production in a purely hadronic scenario 24. l25l l26l. \27\. 



The initial model calculations for J / if) at RHIC |4j yielded interesting predictions for the 
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FIG. 1: (Color online) OPE-based cross sections for J/ip formation and dissociation. 



centrality behavior of the final population Nj/^(tf) normalized by N c5 , or equivalently nor- 
malized by the number of binary collisions. Due to the quadratic behavior of the formation 
rates with N cS , one expects an increase with centrality of formed J/if> per binary collision, 
if the formation reaction dominates. The primary variation in these predictions was due to 
the initial charm quark momentum distributions, which enter through the calculation of the 
formation and dissociation reactivities as defined in Eq. ^ Differences of up to a factor of 
5 or more resulted as the charm quark momentum distributions varied from thermal (pro- 
duced the largest formation efficiency) through perturbative QCD distributions. There may 
be additional uncertainties in overall magnitude. In particular, an implementation of this 
process in a transport model calculation [28j] did not result in a large effect at RHIC energy. 

The first confrontation with experiment involved the initial PHENIX measurements of 
J /if) at 200 GeV Values of dNj/^/dy (y=0) for three centrality bins were extracted, but 
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the most central bin value only yielded an upper bound due to the limitations of statistics. 
Although these limitations prevented a definitive test of the formation model, it was clear 
that the predictions which followed from thermal charm quark distributions were disfavored. 
The pQCD charm quark distributions led to formation rates which were roughly consistent 
with experiment. A parameter scan involving initial temperature, initial charm quark pro- 
duction, transverse expansion velocity of the deconfined region, and initial J/ip production 
number found a surprisingly large region of parameter space within which the formation 
mechanism results were consistent with the PHENIX data points and limits. 

It is interesting that these data are also consistent with binary-scaled J / ip production in 
pp interactions 30] , modified with a small suppression factor. Given the sensitivity of the 
formation process to initial charm production (both magnitude and momentum spectra), 
even the improved statistics anticipated for the PHENIX Au-Au data currently undergoing 
analysis may not be adequate to differentiate between these distinct possibilities. Fortu- 
nately, the momentum spectra for the produced J / ip will also be measured. In anticipation 
of this information, we then proceed to extend the calculations of the kinetic formation pro- 
cess to include the rapidity and transverse momentum spectra of the final J/ip population. 

We return to the expression for in- medium formation of J/ip in Eq.^ modified to calculate 
the momentum spectra of the formed J/ij). Rather than specifying a functional form for the 
charm quark momentum distributions, we will generate a sample of cc events appropriate 
for a given physical situation and then sum over the sample to calculate the differential 
formation reactivity. 



Note that the formation magnitude exhibits the explicit quadratic dependence on N c5 via 
the double summation. 

The formation rate (differential in J/ip momentum) is integrated over time to get the 
spectra of in-medium formed J /if). We assume for simplicity that (a) the spatial density of 
charm and anticharm quarks are equal within the deconfinement region, and (b) the charm 
quark momentum spectra are independent of time. The time-independence assumption 
could in principle be relaxed, at the expense of a more involved calculational scheme. For 
now we simply present the results as limiting cases which bracket a range of possibilities 
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between the two types of charm momentum distributions considered. Then the differential 
J /if) spectra just involve an integral of the inverse deconfinement volume over time, which is 
irrelevant for the normalized momentum spectra. Of course, these spectra will be modified 
by the dissociation process (second term in Eq. [TJ, as well as the spectra of J /if) which are 
initially-produced. We consider the effects of the dissociation process in Sec. II VI 

The differential dependence of the formation cross section which follows from detailed 
balance of the dissociation cross section in Eq. |3] is extracted from the the analog of atomic 
photoionization. In this case the relevant amplitude involves the coupling of a gluon with 
the J /if) color electric dipole of the quarkonium, again in the approximation of a Coulombic 
bound state wave function. The procedure to evaluate the J/if> momentum spectra from 
Eq. 01 just involves a scan for each cc over final state phase space, weighted by the appro- 
priate differential cross section. To determine the sensitivity of the resulting J / if) spectra to 
this particular differential cross section, we have also performed the calculations using an 
isotropic invariant amplitude, normalized to the same OPE-based total cross section. We 
find in general that the effects on the J/if> momentum spectra are negligible. Comparison 
of results will be shown in Sec. IIVI 



III. CHARM QUARK MOMENTUM SPECTRA 
A. Perturbative QCD calculations 

The perturbative QCD spectra are obtained using the results of the full 0{a^) calculation, 



as implemented in the event generator described in |31| . In the case of nuclear beams, the 
parton distribution functions of the proton have been modified to account for nuclear effects, 
usine; the model of ref. The terms of 0(oP s ) corresponding to real-emission diagrams 

will generate a non-trivial transverse momentum distribution for the cc pair. We generated 
samples of unweighted events, which will be used in the following analyses. To ensure that 
all events in the samples have positive weight, we mer ged the negative- weight events arising 



from the cancellation of soft and collinear singularities 31] with the sample of positive weight 
events having a transverse momentum of the cc pair pt{cc) < 0.4 GeV, and smeared their 
Pt(cc) over the — 0.4 GeV range, cc pairs selected from the same event will be labeled as 
"diagonal" pairs: they dominate the J / if) formation in pp and pA collisions, where production 
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of multiple pairs is negligible. In the case of AA collisions, multiple hard scatterings can give 
rise to several cc pairs, and charm (anticharm) quarks from independent pairs are allowed to 
form a J/i(j. These uncorrelated pairs are randomly selected from our sample of unweighted 
events, reproducing in average the correct kinematical distributions of such pairs. 
Shown in Fig. |2] are the shapes of the transverse momentum distributions |q"i|. 
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FIG. 2: (Color online) Transverse momentum distributions from the NLO pQCD cc events at 
RHIC200 energy. The individual charm quark distribution is shown by the circles. The cc pair 
distribution is shown by the squares, and the subset of diagonal cc pairs is shown by the diamonds. 
(Since we use a sample of some 80K cc events, the distribution of all pairs is dominated by the 
off-diagonal pairs.) Also shown by the open triangles is the individual charm quark distribution 
generated using parton distribution functions modified for nuclear effects by the EKS98 |32| factors. 



We plot the spectra of individual charm quarks (with and without inclusion of nuclear 
effects in the PDFs), as well as the transverse momenta of cc pairs. Here and in the following 
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the generic term "cc pairs" will refer to summing over all possible pairings of c and c quarks 
in an event, in the limit of large N c5 multiplicity. The term "diagonal pairs" will refer instead 
to cc pairs coming from the same hard collision. 

The nuclear effects are essentially negligible in the individual charm quark distributions 
(compare the solid circles and open triangles - the scatter at large p t reflects statistical fluctu- 
ations due to our finite sample size). The cc pair distribution is approximately twice as wide 
as that for individual quarks (as measured by (pt 2 ))- This is a result of the independence of 
the transverse momentum vectors (both in size and in azimuthal direction) for off-diagonal 
pairs of c and c quarks produced in different hard scatterings. For diagonal pairs the relative 
azimuthal angle is 7r for all LO events. The small nonzero (pr 2 ) for diagonal pairs thus arises 
entirely from NLO effects. The mean values (pr 2 ) of these distributions are 2.45 GeV 2 for 
the individual quarks, 4.88 GeV 2 for all pairs, and 1.20 GeV 2 for the diagonal pairs. 

Shown in Fig. |3] are the corresponding rapidity distributions. Again the nuclear effects 
are negligible in the individual charm quark distributions. The cc pair distributions are both 
narrower than that for individual quarks. The diagonal pair distribution is somewhat wider 
that that of all pairs, which again follows from the kinematics of diagonal pairs which have 
p t (pair) « from the LO events. These widths can be quantified by the Ay = a/ (y 2 ) — (y) 2 
values, 1.46 for individual charm quarks, 1.25 for diagonal pairs, and 1.03 for all pairs. 



B. Intrinsic kx effects 

We must now modify the charm quark transverse momentum spectra to include effects 
of partonic confinement. This could be accomplished by introducing partonic transverse 
momentum dependence directly into the pQCD calculation. For simplicity, we choose to 
modify the NLO charm quark results by adding an additional transverse momentum 'kick' 
kr to each quark in a cc pair j^. The magnitude is determined by a Gaussian distribution 
with width parameter (kx 2 ), and uniform azimuthal distribution. We will determine the 
width parameter directly from data. For nuclear collisions, the additional effects of initial 
state interactions of the nucleons will be included in this parameterization. 

The resulting charm quark transverse momentum spectra are shown in Fig. 0| for a range 
of (kr 2 )- The width of these spectra satisfy a simple pattern. 

(PT 2 ) ((k T 2 )) = (PT 2 ) ((k T 2 ) = 0) + (k T 2 ) (5) 

10 




FIG. 3: (Color online) Rapidity Distributions from the NLO pQCD cc events at RHIC200 energy. 
Legends are the same as in Fig. |2 

This pattern is again a result of the uncorrelated relative azimuthal angles between the 
transverse momentum kick and the transverse momentum of the original pQCD-generated 
quarks. 

Fig. |S] shows the rapidity distributions of the charm quark events after the (kT 2 ) modifica- 
tion. Aside from small corrections near midrapidity, these spectra are effectively independent 
of (k T 2 ), because the charm quark longitudinal momenta are unchanged by addition of the 
transverse 'kick'. 

We can then investigate the systematic behavior of cc pair distributions as a function of 
(kT 2 )- Fig. |H1 shows the transverse momentum distributions of all pairs and also diagonal 
pairs. The symbol legends are the same as that in Fig. 0] and Fig. |S] The distributions for 
all pairs (dominated by the off-diagonal pairs) shown by the corresponding solid curves are 
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FIG. 4: (Color online) Charm quark transverse momentum distributions including intrinsic 

significantly broader than those for the diagonal pairs at (k T 2 ) = because the pair px = 
for the LO events. However, that difference decreases as (/ct 2 ) increases, because the quarks 
in the off-diagonal pairs receive uncorrelated fc^-kicks. This is verified by a linear fit to the 
(Pt 2 ) of the pair distributions. For the full range of generated diagonal pairs, one obtains 



to be compared with 



(Pt) diagonal cc 



(pT 2 )all cc 



1.2GeV 2 + 4(fc T 2 ), 



4.9 GeV 2 + 2(k T 2 ) 



(6) 



(7) 



for the sample with all pairs. 

It is also evident that the diagonal pair sample exhibits some scatter at the larger px 
values due to statistical fluctuations limited by the total number of generated events which 
is of order 10 4 . This limitation does not appear until much higher p? for the all pairs sample 
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FIG. 5: (Color online) Charm quark rapidity distributions including intrinsic fey. 

with total number of order 10 8 . The effect is evident if we truncate the high-pr events in 
the diagonal pair sample. The resulting width satisfies 



(pT 2 )dia 9 onai Apt < 5GeV) w 0.8 GeV 2 + A(k T 2 ). (8) 

This parameterization will be useful in comparison with initial PHENIX data in pp and 
d-Au interactions, where the measured pt of J/ ip are limited to the pt range in Eq. |HJ 

The rapidity spectra of the cc pairs are shown in Fig. again using the same symbol 
legends to indicate (^t 2 )- One sees that there is not a dramatic effect in the range of (k T 2 ) 
up to 2.0 GeV 2 . As (kx 2 ) increases, the width Ay decreases by less than 10 % for both all 
pairs and diagonal pairs. 
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FIG. 6: (Color online) cc pair transverse momentum distributions including intrinsic kx- Shown 
are all pairs (solid lines) and diagonal pairs (symbols, as in Fig. [SJ. 

C. Thermal plus flow charm quark distributions 

We also have generated a set of charm quark momenta to investigate the opposite scenario, 
in which the interaction with the partonic constituents in a deconfmed region is sufficiently 
strong to result in their thermalization in the expanding medium. This scenario is not favored 
by conventional arguments based on small heavy-quark cross sections and the inhibition of 
medium-induced radiation due to the dead-cone effect . Analysis of high p T leptons in 
Au-Au collisions as a measurement of charm quark production are consistent with spectra 
unmodified from the initial pQCD calculations |3(|. However, it was subsequently shown 
[37 | that the resulting spectra of heavy-quark mesons and electrons for < 3 GeV cannot 
distinguish between the two different scenarios. Recently the possibility of heavy-quark 
thermalization was found to be accelerated by interactions involving resonant hadronic states 
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FIG. 7: (Color online) cc pair rapidity distributions including intrinsic kj>. Shown are all pairs 
(solid lines) and diagonal pairs (symbols, as in Fig. EJ. 

in a plasma 38]. In any event, it is interesting to consider such a heavy-quark momentum 
spectrum in the kinetic formation process to illustrate the impact of this very different 
scenario on the resulting bound state spectra. 

To generate this set of charm quark momenta, we use the distribution which follows from 
numerical results of transverse-boosted Bjorken boost invariant hydrodynamics 



dN 



" R j t ( PTsmhy T (r)\ (m T coshy T (r) 
oc m T I rdrlo [ — ) Ai 1 



dpT 2 ^"^ Jo T J" L \ T 7' (9) 

where the transverse mass rrix = \/m 2 + pr 2 and R is the transverse system size. The 
hadronization temperature T and the rapidity of transverse expansion yr{r) are deter- 



23 determined T = 128 MeV and 



mined from a fit to the light hadrons. Ref. 
yr(r) = tanh -1 ^^?")) with a linear boost profile for the transverse expansion velocity, 
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Pr{f) = Pmax(ji) with j3 max = 0.65 for central collisions. This form has also been used to fit 
the transverse momentum spectrum of hadrons in a statistical model in which the freeze- 
out and hadronization temperatures are identical 0|. In this case the linear profile of the 
transverse expansion was taken to vary with transverse rapidity rather than velocity, but for 
typical fimax values the difference is not significant. For the calculations here, we have used 
the linear profile in transverse rapidity. We also consider an alternate set of parameters, 
T = 170 MeV and (3 max = 0.5 which emerge from a coalescence study of thermal partons 
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which also fit the light hadron spectra at RHIC 

Fig. |H1 shows the transverse momentum spectra of charm quarks which result from Eq. El 
using the parameter sets (T, y Tmax ) = (128 MeV, 0), (170 MeV, 0), (128 MeV, 0.65) and (170 
MeV, 0.50). For no flow, the distributions for different temperatures are of course different. 
However when augmented by the corresponding transverse flow the resulting distributions 
are almost identical, since both yTmax values have been adjusted to fit the observed flow at 
the corresponding freezeout temperatures. 

The pQCD distributions shown for comparison are much broader and have an obvious 
non-thermal behavior for large pp- The widths can be characterized by the calculated (pt 2 ) 
values of 0.47, 0.67, 1.35, and 1.24 GeV 2 for the thermal and thermal plus flow distributions, 
compared with 2.45 and 3.45 GeV 2 for the pQCD examples. 

The transverse momentum spectra of the thermal plus flow cc pairs are shown in Fig. El In 
this case there is no set of diagonal pairs since all identity has been lost in the thermalization 
process. The pair distribution widths are twice those of the corresponding distributions for 
single charm quarks, as expected due to the uncorrelated relative azimuthal angles. Their 
{Pt 2 ) values are 0.93, 1.32, 2.70, and 2.46 GeV 2 , respectively. Again these are much smaller 
than the pQCD comparison distributions, which have (pt 2 ) of 4.9 and 6.9 GeV 2 . 



IV. DETAILS OF KINETIC FORMATION MODEL CALCULATIONS 

We generate the transverse momentum and rapidity spectra of the formed J/ ip according 
to Eq. HI For each cc pair, we transform to the pair CM-system and weight the formation 
event by v re \ op. The formation cross section op is obtained via detailed balance from the 
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FIG. 8: (Color online) Thermal plus flow charm quark transverse momentum distributions. Shown 
are thermal without flow (solid and dashed lines), thermal with flow (star and dotted line), and 
pQCD (symbols). 



dissociation cross section ajj, 



= 48 (g-M J/v , 2 ) 2 
36 s (s — 4m c 2 ) 



(10) 



in terms of the charm quark mass m c and the energy invariant s. <jp> is taken from Eq. 01 
with modified threshold gluon momentum to account for the finite mass of the J/ip. 

The angular distribution of the J / ip in the cc rest frame is taken from the corresponding 
behavior of atomic or nuclear photodissociation which gives 



d 3 p 



oc sin 2 9. 



(11) 



The sensitivity to the high energy behavior of ap is probed by imposing a maximum energy 
cutoff parameter d max |52| which limits the invariant mass of the cc pairs which can initiate 
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FIG. 9: (Color online) Thermal plus flow pair transverse momentum distributions. Shown are 
thermal without flow (solid and dashed lines), thermal with flow (star and dotted line), and pQCD 
(symbols). 

the formation of a J/ip: 

V 1 Sec < Mj/f + d max (12) 

The resulting rapidity and transverse momentum spectra are given by the black circles 
compared in Fig. (TUl and Fig. [TTJ respectively. 

Our choice of initial default parameters set the binding energy eo = 0.63 GeV, the mass 
cutoff parameter d max = 3.0 GeV, and no kick ((hr 2 ) = 0). The curves refer to a sample 
of 10 4 independent cc pairs, resulting in 10 8 possible pairings. 

The impact on these distributions of the various dynamical assumptions is shown by three 
additional curves: 

• J/i/j, isotropic matrix element (empty squares): the angular dependence in Eq. fTTI is 
taken to be flat; 
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FIG. 10: (Color online) Comparison of rapidity spectra for cc pairs and formed J /if), all using 
default parameters. 



• weighted cc pairs (solid lines): we plot the kinematical variables of the cc pairs weighted 
by the J/ip formation probability v re i o>, which differs from the formation process by 
elimination of exact kinematics relating a cc pair to the spectrum of J / ip which follow 
from the in-medium formation reaction; 

• cc pairs (dashed lines): we plot the kinematical variables of the cc pairs which result 
from all combinations of charm quarks as produced in the pQCD processes. 

One sees that for both spectra the J / ip formation is not sensitive to the angular dependence 
of the amplitude, as evidenced by the equality of the circle and square plots (except a small 
effect at the largest values of transverse momentum). The primary effect of the formation 
dynamics is to produce a narrowing of all spectra, compared with that of the cc pairs (dashed 
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FIG. 11: (Color online) Comparison of transverse momentum spectra for cc pairs and formed J/ip, 
all using default parameters. 

line). For the rapidity spectra, even the weighted cc pairs (solid line) are identical to the 
Jji>. 

The effects of intrinsic (/ct 2 ) on the J/ip spectra are shown in Fig. IT2l and Fig. ITBl 

The rapidity spectra are essentially independent of (&t 2 )- The transverse momentum 

spectra increase with (kx 2 ), as expected. We can parameterize this effect in terms of (pt 2 ) 

for the formed J/ip, as 

(p T 2 )j/<p =2AGeV 2 + (k T 2 ). (13) 

It is interesting to compare this form with the corresponding behavior of all cc pairs in 
Eq. [7J in which the value at (kx 2 ) = is larger (4.9 GeV 2 ) as is the rate of increase (2). 
A fit of an intermediate case, using all cc pairs weighted with the formation cross section 
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FIG. 12: (Color online) Effects of intrinsic {kx 2 ) on the J/ip formation rapidity spectrum. 
Vrei &f yields an intermediate result, 



(PT 2 ) ( 



3.2 GeV 2 + 1.5 (k T 2 ). 



(14) 



/ ccbar (weighted) 

Next we show the sensitivity of the formation results to variations in some of the default 
parameters. Fig. El shows the results for the normalized transverse momentum spectra. 
The first three curves show variation in d max (symbols). The final two curves (solid and 
dashed lines) show the effects of reducing the binding energy to eo = 300 MeV, more than 
a factor of two below the vacuum value. In addition, the last curve has the quark mass 
increased to m c = 1.7 GeV, which restores the J/ip mass to its vacuum value. One sees 
that overall these changes produce negligible variations in the pt spectra. Lastly, we look 
at extreme variations in the functional form of the cross sections to determine the degree of 
stability of the J/ip formation p T spectra. Fig. ITBI shows some of these effects. The effect of 
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spectrum. 

omitting an energy cutoff (d max — > oo) is shown by the thin solid line, which is essentially 
unchanged from the default parameter result (shaded circles) with d max = 3 GeV. This 
might be understood from the behavior of the OPE-cross sections, since the entire region 
around the peak values is generally contained within the default parameter choices. We next 
consider the effect of a cross section independent of y/s, shown by the dashed line. (The 
calculation is virtually unchanged whether one chooses <jp or <Jd to be constant). This choice 
of cross section produces only a small change in the J /if} transverse momentum spectrum. 
The most radical departure is to use constant cross sections and no energy cutoff. These 
results are shown by the solid line and circle, and do produce a significant change in the 
Pt spectrum. We would certainly not expect, however, that the J ftp formation probability 
would remain independent of the cc pair relative momentum throughout the entire range 
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FIG. 14: (Color online) Effects of parameter variations on the J/ip formation transverse momentum 
spectrum. 



populated by the pQCD initial production process. To the contrary, one would expect that 
the large invariant-mass pairs would preferentially hadronize into the higher-mass states in 
the open charm spectrum. 

Finally, we investigate the effects of in- medium dissociation on the formed J ftp spectrum. 
This dissociation will influence the final pt spectrum in two competing scenarios. In the 
first, one expects that high-p^ J/ip which are initially produced will originate preferentially 
in central collisions, where the Cronin effect is most efficient, and thus be subject to a large 
suppression. However, one also expects that the J / ip with large px will preferentially escape 
the deconfinement region. In a recent study [46] incorporating these effects for J ftp propa- 
gating in a region of deconfinement at both SPS and RHIC energy, it was found that this 
second "leakage" effect dominates. We will use numerical values taken from initial studies of 
this effect J43 1 and updated versions 44, 45] in the form of pr-dependent suppression factors 



23 



10 




I I I I I I I I I I I 

2 4 6 8 10 

P T (GeV) 

FIG. 15: (Color online) Dependence of the J/ip formation transverse momentum spectrum on the 
formation cross section. 

for J/ip initially present in an equilibrating and expanding parton gas. For the application 
to our formation spectrum, we need a suppression factor for J/ip formed at a continuum of 
initial times. However, the time of formation is biased toward early times when the charm 
quark density is maximum. Hence we use two typical suppression factors from Ref. j^ij], 
which have their minimum value at px = and approach unity for pt ~ 10 GeV. These 
are shown in Fig. The effect of these suppression factors on the formation spectra is 
shown in Fig. El One sees that the normalized formation spectra are suppressed in the 
low-pr region by amounts from 10 - 25%. However, the normalized spectrum including both 
formation and suppression is virtually unchanged in the low-pr region (which dominates the 
normalization factor), and the entire effect of dissociation is pushed out to very large pt- 
Hence the normalized formation spectra alone provide sufficiently accurate predictions, at 
least in the region < px < 6 GeV where the majority of J/ip are found. This effect can 
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FIG. 16: (Color online) In-medium J ftp suppression factors. 



be quantified by the small increase in (pr 2 ), from 3.59 GeV 2 to 3.69 GeV 2 



V. PREDICTIONS FOR AU-AU INTERACTIONS AT RHIC 



A. Charm quark distributions from pQCD 

In order to make predictions of the kinetic formation model for J ftp at RHIC, we need 
to fix the {kr 2 } parameter. We start with J ftp production in pp interactions at 200 GeV. 
There are PHENIX data on the rapidity and transverse momentum spectra from both run2 
I30I ] and run3 j^], which we exhibit in Fig. [THl Since multiple cc pair production in pp 
interactions is negligible, we consider only the diagonal cc pairs. These are shown by the 
solid curve, and are seen to describe the data reasonably well. (Recall from Fig. that 
the cc pair rapidity spectra are essentially independent of (kr 2 )-) One can interpret this 
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FIG. 17: (Color online) The in-medium J ftp normalized formation plus suppression spectrum 
shown by the small filled squares is compared with the normalized formation-only spectrum shown 
by the large open squares. 

agreement as consistency with any mechanism in which J/ip is produced via hadronization 
of initially-produced cc pairs. Also shown by the dotted curve are the pairs modified by the 
in-medium formation probability, proportional to ap. This process of course cannot occur 
in a pp interaction, and it is gratifying to see that it is incompatible with the data. 

The transverse momentum spectrum is shown in Fig. EH The PHENIX dimuon data are 
presented for both positive and negative rapidity, which have smaller uncertainties than the 
central dielectron data. The solid triangles are the pQCD cc diagonal pairs, which of course 
must be augmented by some (kx 2 ) to include initial state and confinement effects. The solid 
line and adjacent circles and diamonds represent (h? 2 ) = 0.5 ± 0.1 GeV 2 , and provide an 
acceptable fit to the data (fluctuations at large pt are statistics limited from the number of 
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FIG. 18: (Color online) Rapidity spectra of J/ip production in pp interactions at 200 GeV. 

generated diagonal pQCD events). 

There are also preliminary data on J / ip production in d-Au interactions (^tJ . The width 
of the pt spectra is found to be larger than that measured in pp interactions. 



(p T 2 )d-Au 



(PT 2 ) PP 



(15) 



1.77 ±0.35 GeV 2 {y = -1.7) 
1.29 ± 0.35 GeV 2 (y = +1.8) 
This broadening of pt distributions for particles produced on nuclear targets is well 
known, and fits a natural pattern which emerges from initial state elastic scattering of a 
projectile in the nuclear target For p-A (or d-A) interactions, the increase can be 

expressed clS db change in (pt 2 ), as 



{Pt 2 ) p a ~ (Pt 2 } 



pp 



A [n A -l], 



(16) 
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FIG. 19: (Color online) Transverse momentum spectra of J /ip production in pp interactions at 200 
GeV. 

where ha is the impact-averaged number of inelastic interactions of the projectile in nucleus 
A, and A 2 is the square of the transverse momentum transfer per collision. For a nucleus- 
nucleus collision, the corresponding relation is 

(pt 2 )ab - (pt 2 ) pp = A 2 [n A +n B - 2]. (17) 

We use the measured J/ip broadening in d-Au to determine the appropriate (kx 2 ) value for 
the pt distribution of initially-produced diagonal cc pairs through Eq. |U]or|Hl with the result 
(kT 2 )d-Au — (k T 2 ) pp = 0.4 ± 0.1. (Since the measured values for d-Au interactions differ 
between positive and negative rapidity, we use their average to partially compensate for the 
existence of final-state effects.) Finally, using (/ct 2 ) pp and (kT 2 )d-Au extracted from data, 
in combination with Eq. IT7I leads to (kr 2 ) Au-Au = 1.3 ± 0.3 GeV 2 . We use this value to 
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determine the pt distribution of initially-produced charm quarks in Au-Au interactions. It 
is interesting to also look at the equivalent parameters which follow from Eq. and Eq. 
[TBI One finds = 5.4 for minimum bias d-Au interactions at RHIC energy (using a pp = 
42 mb), which leads to A 2 = 0.35 ± O.lAGeV 2 . We note that the relatively large uncertainty 
comes entirely from the difference in p? broadening in d-Au between positive and negative 
rapidity. It is interesting to see the energy dependence by comparing this value with that 
extracted from J/ip pA data at fixed-target energy 48(, A 2 = 0.12 ± 0.02GeV^ 2 . 

Shown in Fig. [201 are the predicted rapidity spectra of J / ip in Au-Au interactions at 200 
GeV. The solid triangles use initially-produced diagonal cc pairs with (fcr 2 ) = 1-3 GeV 2 , 
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FIG. 20: (Color online) Predicted rapidity spectra of J ftp in Au-Au interactions at 200 GeV. 
Triangles are for initial production via diagonal cc pairs. Circles are for in-medium formation via 
all pairs. Sensitivity of the formation spectrum to variation of (kx 2 ) within the range (&t 2 ) 
1 ..'! ± 0.3 GeV 2 is indicated by the spread in the solid lines. The corresponding spread in the lines 
for diagonal pairs covers the range < (A>r 2 ) < 2 GeV 2 . 
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and the various lines indicate the small sensitivity to variations of (&t 2 ) Au-Au within un- 
certainties. This distribution is also virtually the same as those using (kT 2 ) pp , which fit 
the pp data as shown in Fig. EI The solid circles and associated lines show the results 
of a formation calculation using all cc pairs which can recombine in the medium. These 
spectra are substantially narrower, and provide a prediction which signals the existence of 
the formation mechanism. Since the formation process is largest for central collisions, where 
dissociation of the initial production yield should be most efficient, one would expect that 
the total J/ip rapidity spectrum should change from narrow for the most central collisions 
to wide for the most peripheral collisions. 

Fig. 1211 shows the corresponding transverse momentum spectra. The solid triangles again 
use initially-produced diagonal cc pairs with (kx 2 ) Au-Au = 1-3 GeV 2 . Note that this dis- 
tribution is wider than those using {kT 2 ) pp , which fit the pp data shown in Fig. The 
solid circles and associated lines (which include the uncertainty in {kx 2 ) Au-Au) show the 
results of a formation calculation using all cc pairs which can recombine in the medium. 
These spectra are substantially narrower, and provide another prediction which signals the 
existence of the formation mechanism. The same consideration of centrality dependence as 
presented for the rapidity spectra also apply to these transverse momentum spectra. 

B. Charm quark distributions from thermal plus flow 

The calculation of J/ip follows from the same cross section which has been used for the 
pQCD charm quark distributions. We use generated cc events as calculated in Sec. IIII C\ 
where it was demonstrated that the charm quark transverse momentum distributions which 
follow from thermal plus flow parameters T = 128 MeV, yrmax = 0.65 and T = 170 MeV, 
VTmax = 0.50 are essentially identical. 

We show in Fig.|221the transverse momentum spectra evolution which follows from charm 
quarks with T = 170 MeV, i/Tmax = 0.50. One sees a downward progression in widths, 
starting with (pr 2 ) = 2.5 GeV 2 for cc pairs, to 2.0 GeV 2 for weighted cc pairs and 1.7 GeV 2 
for J I ip formation. It is important to note that the concept of diagonal pairs no longer exists 
for this set of charm quarks, since their formation identity has been entirely erased during 
the thermalization process. As might have been expected with a thermal distribution, these 
widths are all smaller than any of those which follow from pQCD charm quark distributions. 
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FIG. 21: (Color online) Predicted transverse momentum spectra of J/ip in Au-Au interactions 
at 200 GeV. Triangles are for initial production via diagonal cc pairs. Circles are for in-medium 
formation via all pairs. Sensitivity of the formation spectrum to variation of (£>r 2 ) within the range 
(kx 2 ) = 1.3 ± 0.3 GeV 2 is indicated by the spread in the solid lines. 

This is shown in Fig. |2H1 The filled diamonds and triangles are the J/ip formation and 
diagonal cc pair distributions which follow from the pQCD quark distributions as calculated 
in Sec. IV Al The solid squares and dashed line are for the two equivalent thermal plus 
flow charm quark distributions, and for completeness a formation spectra without charm 
quark transverse flow is shown by the solid circles. It is clear that in-medium formation of 
J/i/j starting from either the pQCD or thermal plus flow charm quark distributions predict 
transverse momentum widths (characterized by (pr 2 )) which are markedly smaller than that 
for initially-produced J/ip. The difference between the two formation possibilities would be 
reflected by the significant difference between the shapes of the pr distributions. 
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FIG. 22: (Color online) Transverse momentum spectra comparison for J/tp and ccbar pairs for a 
thermal plus flow initial quark distribution. 

In Fig. |2Uwe compare these results with alternate model calculations in the literature. 
The filled diamonds and squares are the J/ip transverse momentum spectra from our in- 
medium formation calculations using the two extreme assumptions about the charm quark 
momentum distributions, pQCD (diamonds) and thermal plus flow (squares). We compare 
this with a calculation which uses a coalescence model to produce J lib from charm quark 

n 

pairs |49[. The solid line is the result using charm spectra generated by Pythia (compare 
with our pQCD formation), and the dashed line uses thermal plus flow charm spectra in 
the coalescence model (compare also with our formation using thermal plus flow). One sees 
that the systematics of the model spectra are the same, but the coalescence process yields 
somewhat broader spectra than our direct in-medium formation. Also shown by the filled 
circles is the J/ib spectrum which follows from a statistical model of hadronization. Here 
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FIG. 23: (Color online) J/ip formation transverse momentum spectra dependence on initial charm 
quark distributions. 

the J/ip would inherit a thermal plus flow kinetic distribution This distribution is 

broader than the underlying charm quark thermal plus flow case, but narrower than the 
pQCD formation and coalescence processes. The width of the J/tp transverse momentum 
spectra at RHIC has also been predicted in j^. The value quoted for the most central 
collisions is (pr 2 )= 3.3 GeV 2 for the most central collisions. For comparison, our pQCD 
formation process predicts 3.6 GeV 2 in a region where it is the dominant process. 



VI. SUMMARY 



We have shown that in-medium formation of heavy quarkonium states utilizing charm 
and anticharm quarks coming from independent hard scatterings results in normalized mo- 
mentum spectra which retain a memory of the underlying quark distributions. In the case 
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FIG. 24: (Color online) J ftp formation transverse momentum spectra comparison with alternative 
models. 

that the quark distributions follow from initial production in pQCD processes, both the 
rapidity and transverse momentum spectra of the formed heavy quarkonium will be nar- 
rower than those expected from diagonal pairs in the absence of a color-deconfined medium. 
Explicit calculations for J / ip formation in Au-Au collisions at 200 GeV are performed, using 
initial data in pp and d-Au interactions to fix some parameters. A striking feature is the 
non-monotonic behavior of the transverse momentum spectrum widths (pr 2 ) in the progres- 
sion pp to pA to AA. In the absence of in-medium formation, one would expect (pr 2 ) to 
increase monatonically with the colliding system size, due to the initial state confinement 
and nuclear broadening effects. Our calculations for in-medium formation in AA collisions 
predict (pt 2 ) substantially smaller than expected if the J/ip were due to initial production 
(see Fig.EU) In fact, the predicted {pt 2 )aa = 3.59 GeV 2 is even smaller than the measured 



value in d-Au interactions, (pr 2 ) 



d-Au 



4.25 GeV 2 . 
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This effect is maximum for central AA collisions, and will revert to the "normal" behavior 
for peripheral collisions. Uncorrelated pairs, composed of one quark and one antiquark which 
do not originate from the same initial interaction, have an invariant-mass spectrum which 
becomes harder for larger p? of the pair. Given that the charmonium formation dynamics 
favors lighter invariant masses, we expect that the contribution of uncorrelated pairs will be 
more dominant at smaller pr, leaving a signature in the overall pt slope. The softening of the 
Pt spectrum for uncorrelated pairs, in other words, is the result of an increased contribution 
to charmonium production at lower px- 

We have also considered the formation process using thermal charm quarks which flow 
with the expanding medium. The resulting spectra are substantially narrower and retain a 
form which reflects that of the underlying heavy quarks. Overall, our normalized spectra 
appear to be quite robust with respect to variations of the model parameters, independent 
of the absolute magnitude of in-medium formation. 
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